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Abstract 

A model independent formalism for the electron positron annihilation reaction e + + e~ — > 
ai(1260) + 7r has been derived. The differential and total cross sections and the elements of the 
spin-density matrix of the a\ -meson were calculated in terms of the electromagnetic form factors of 
the corresponding 7*ai7r current. Simple models of a\ form factors have been fitted to the available 
cross section data and they allow to give numerical predictions for the different observables. 

PACS numbers: 12.20.-m, 13.40.-f, 13.60.-Hb, 13.88.+e 
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I. INTRODUCTION 



The electron positron annihilation into hadrons constitutes an important source of infor- 
mation on the internal structure of the mesons: the light quarks and their interactions as 
well as the spectroscopy of their bound states. The experimental data about these reactions 
in the low-energy region are also relevant to the determination of the strong interaction 
contribution to the anomalous magnetic moment of the muon, to the test of standard model 
predictions for the hadronic tau-lepton decay, which is related by the conservation of vector 
currents. 

Recently, the construction of the new detectors with a large solid angle, which can operate 
at new colliders with high luminosity opened new possibilities for the investigation of the 
reactions e + + e~ — > multihadrons Not only the statistic is highly increased, but also 
the possibility to detect charged as well as neutral pions allows to draw conclusions on the 
nature of the intermediate states. 

In the energy region 1 < W < 2.5 GeV (W is the total energy of the colliding beams) 
the process of four pion production is one of the dominant processes of the reaction e + + 
e~ ^hadrons. Its cross section is larger than 2ir production and comparable to e + + e~ — > 
fi + + uT. 

The process of e + e annihilation into four pions was firstly detected in Frascati [2] and 
later on in Novosibirsk [3j. Through a simultaneous analysis of the differential distributions 
in two final channels: 27r + 27r _ and 7r + 7r~27r°, it was shown in [l] that the reaction pre- 
dominantly occurs through the ai(1260)7r and ujtc intermediate states in the energy range 
1.05-1.38 GeV. It was also found that the relative fraction of the ai(1260)7r state increases 
with the beam energy. The measurement of the e + e~ — > 7r + 7r~7r + 7r~ cross section was ex- 
pended to lower energies. Data obtained with larger statistical and systematic precision 
4j confirmed that the dominant production mechanism is consistent with the a!(1260)7r 
intermediate state. 

The process of the multihadron production at large energies was also investigated with the 
BABAR detector at the PEP-II asymmetric electron-positron storage ring using the initial- 
state radiation [5|]. In particular, the cross section for the process e + e~ — > n + n~n + n~ was 
measured for center-of-mass (CMS) energies from 0.6 to 4.5 GeV, providing evidence of 
a resonant structure, with preferred quasi-two-body production of ai(1260)7r. A detailed 
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understanding of the four-pion final state requires also information from final states such as 
7r + 7T - 7r°7r , to which the p + p~ intermediate state can contribute. A summary of the hadronic 
cross section measurements performed with B ABAR via radiative return is given in Ref . {(| . 
In this paper we consider the reaction 

e-(A: 1 )+e + (A; 2 )^ai(p 1 )+7r(p 2 ), (1) 

where a± is the axial-vector meson ai(1260) with the following quantum numbers I G (J PC ) = 
1~(1 ++ ). The notation for the particle four-momenta is given in brackets. 

The determination of the pseudoscalar-meson FFs as pions and kaons requires only cross 
section measurements. They have been extensively studied both in the space-like and time- 
like regions (see, for example, [7]). Light vector mesons are less known, because their 
experimental determination is more difficult, due to their short lifetimes. However, the t 
dependence of the cross section for diffractive vector-meson electroproduction gives (model- 
dependent) information on the charge radius, and radiative decays, such as p + — > 7r + 7r°7, 
allow to obtain their magnetic moment. 

From the theoretical point of view, the processes of the vector- and axial-mesons pro- 
duction in the electron-positron annihilation were considered in a number of papers. The 
predictions for the differential and total cross sections of the reaction e + e~ — > af + n T was 
given in Ref. 8j in the framework of the hard-pion current algebra models. 

It was shown that the reaction cross section alone could be, in principle, discriminative 
toward models. Using VMD model, the authors of Ref. [s| investigated the reaction e + e~ — > 
mesons assuming two-body (or quasi-two-body) final states, as ai(1260)7r and p + p~. All 
FFs were taken equal to unity. Estimations of the cross sections of the processes e + e _ — > 3n, 
47r were also obtained, using the VMD model, in Ref. [lo| . 

It appears that the magnitude and energy dependence of the cross section alone can not 



constitute a decisive test on the validity o: 
in Ref. 11], where existing models [12 



' VMD models. A recent discussion can be found 
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14| have been phenomenologically modified, 



including parameters to be fitted on the data. A good description of the cross section is 
obtained assuming aiir intermediate state, in addition to p and 7r, and including higher p 
resonances for energies over 1 GeV. 

Due to the conservation of vector current the cross section of the e + e~ — > An process can 
be related to the probability of the r — > Ai\v T decay. Therefore, all realistic models describing 



3 



;he first process, should also be applicable to the description of the latter one. It was found 
I5I that the assumption of the ai(1260)7r dominance is in qualitative agreement with all 
available data. The free-parameter investigation of the branching ratios and distribution 
functions of the four particle decay of r — > prniu, in terms of the effective chiral theory of 
mesons, is consistent with the data [if]]. The theory predicted the ax dominance in these 
four particle decay of the tau-lepton. 

The purpose of this paper is to calculate the differential (and total) cross sections and the 
elements of the spin-density matrix of the ai-meson in terms of the electromagnetic form 
factors (FFs) of the corresponding 7*ai7T current. A model independent formalism, derived 
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allows 



171 ] for spin one particles, and applied to the process e + + e~ — > p + + p~ in 
to express the experimental observables (differential cross section, polarization observables, 
elements of the density matrix..) in terms of hadron FFs. In annihilation reactions, these 
FFs should be known, or extrapolated from the space-like region into the time-like (TL) 
region, on the basis of analytical arguments. 

II. FORMALISM 

The following derivation is based on the one-photon exchange mechanism. In principle, 
at large q 2 , one should take into account the two-photon-exchange contribution, as it was 



suggested a few decades ago jl9|] . In case of spin one particles, model independent properties 



of the two-photon-exchange contribution in elastic electron-deuteron scattering have been 



discussed in Refs. 
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2 if ] . However no experimental evidence has been found, up to now, 



on the presence of two-photon-exchange in the scattering 22| as well as in the annihilation 



channels 



231 ]. and we do not include such contribution which makes the formalism very 



complicated. 



A. The spin structure of the matrix element 



In the one-photon approximation, the differential cross section of the reaction (OQ) in terms 
of the hadronic, W^ u , and leptonic, L^ u , tensors, neglecting the electron mass, is written as 

da a 2 p 



dVL q e 2W 



L^W^, (2) 
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where a = 1/137 is the electromagnetic constant, p = y(q 2 + m 2 — M 2 ) 2 — Am 2 q 2 /2W is 
the final-particle momentum in the reaction CMS, m and M are the masses of the pion and of 
the cxi meson, respectively. The four momentum of the virtual photon is q = ki + k 2 = P1+P2, 
with q 2 = W 2 , and W is the total energy of the initial beams (note that the cross section is 
not averaged over the spins of the initial beams). 

The leptonic tensor (for the case of longitudinally polarized electron beam) is 

L(i V = ~q 2 g^u + 2(k lli k 2v + k 2fl k lv ) + 2i\e fiuap k la k 2 p , (3) 

where A is the degree of the electron beam polarization (further we assume that the electron 
beam is completely polarized and consequently A = 1). 

The hadronic tensor can be expressed via the electromagnetic current J M , describing the 
transitions 7* — > irai, as follows 

w, u = (4) 

The hadron tensor W pv can be expressed in terms of FFs of the 7* — > ircii transition, 
using the explicit form of the electromagnetic current J M . The spin-density matrix of the 
a%— meson is composed of three terms, corresponding to unpolarized, vector and tensor 
polarized meson: 

Pfiv — - \ 9tiv -^p~ J + -^-^iivpoSpPla + O^pu- 

Here and are the ai— meson polarization four vector and quadrupole tensor, re- 
spectively. The four vector of the a\— meson vector polarization and the a\— meson 
quadrupole-polarization tensor Qp U satisfy the following conditions: 

s 2 = -1, spi = 0, Qf, u = Q Ufl , Qpp = 0, pipQp V = . 

Taking into account Eqs. (TjJ and the hadronic tensor in the general case can be 
written as the sum of three terms 

= VM°) + W^(V) + W^(T), (6) 

where Wp V (0) corresponds to the case of unpolarized particles in the final state and 
Wp U (V)(Wp U (T)) corresponds to the case of the vector (tensor) polarized a x — meson. 

These expressions are general, for any spin one particle in the final state. Let us consider 
more particularly, the reaction e + + e~ — > n + a\ which has been shown to be the main 
contribution to the 4it final state. 
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The electromagnetic current of the 7* — > 7rai transition is described by two FFs. Assum- 
ing the P and C-invariance of the hadron electromagnetic interaction this current can be 
written as IOJ 



J, = h(q 2 )(q 2 u; - q ■ u* q ,) + f 2 ( q 2 )( q ■ P2 u; - q ■ u* P2fi ), (7) 

where U p is the polarization four- vector describing the spin one ai-meson, and fi(q 2 ), (i = 
1, 2) are the electromagnetic FFs describing the 7* — > irai transition (note that we singled 
out explicitly the electron charge e from the expression for the electromagnetic current). 
The FFs fi{q 2 ) are complex functions of the variable q 2 in the region of the TL momentum 
transfer (q 2 > 0). 

In case of real photon, fi does not contribute, and the value /2(C)) can be obtained from 
the experimental data on the decay width T(ai — ► Try). The expression of the radiative decay 
width T in the axial-vector meson rest frame is 

" 1 W 'M\ 2 , (8) 



tt24M 2 

where u is the photon momentum in the axial-vector meson rest frame, u = (M 2 -m 2 )/2M, 
and the matrix element of the radiative decay is written as: M. = ee^J^, where e M is the 
photon polarization four-vector. So, the value of FF / 2 at point q 2 = can be obtained 
from the following expression 

FA/f 3 



Another estimation of these FFs was obtained in Ref. [10| where the processes e + e~ — > 3tt, 
Att were studied. The spin structure of the matrix element of the aip-n transition is similar 
to the 7*ai7r transition. The two coupling constants describing the aipn transition were 
estimated by using the method of hard pions and knowing the width of the ai-meson (taken 
about 130 MeV) 10|. The corresponding values: g 2 /Air w 1.3 and (3 = 1.6 were obtained. 
The following relations hold: 

ef 1 (q 2 = m 2 p ) = ^-, ef 2 (q 2 = m 2 p ) = -L(0 - 1). 

The explicit expressions for the contributions to the hadronic tensor are: 
- Unpolarized term 1^,(0): 
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W pu (0) = W 1 (q 2 )g pu + 



W 2 (q 2 ) „ _ _ q p q v _ Pi<? 

-pi/xPii/ , = V — , Pin = Pi» - —% 



M 2 



where 



W^q 2 ) = -\q 2 h + \{q 2 - M 2 + m 2 )f 2 \ 2 , W 2 {q 2 ) = q 2 (q 2 \h + / 2 | 2 ~ M 2 |/ 2 | : 



(10) 



- Term for vector polarization W pv {y): 



W^uiV) = -^V 1 (q 2 )e puap s a q p + ^V 2 (q 2 ){p lp e„ aap s a q a p lp - p lu e paap s a q a p lp ] + 
+ jpV 3 (q 2 )[p lp e uaap s a q CT pi p + P\ v e paap s a q a p Xp \ 
V x {q 2 ) = —M 2 (q 2 + M 2 — m 2 )~ l \q 2 fi + ^(g 2 — M 2 + m 2 )f 2 \ 2 , 



2\-l 



V 2 (q 2 ) = —Mq(q + M — m 

I( 3g 2_ M 2 + m 2 )Re/i/ , 

V 3 (q 2 ) = -l-M'qHmfj;, 



g 2 \fi\ 2 + ^ 2 -M 2 + m 2 )\f 2 \ 2 + 



(11) 



Term for tensor polarization W pu (T): 



W pv (T) = T 1 (q 2 )Qg^ + T 2 (q 2 )-^p lp p lu + T 3 (q 2 )(p lp Q u +p lu Q^ + 



M 

T^q^Q^ + iT b (q 2 )(p lpL Q v - p lv Q 



Hit 



(12) 



where 



Uq 2 ) 



Quuqu nQi Qu,qii 0, Q uv Qui/ A Q r, Qua r, Qvai 

q 2 g 4 q 2 q 2 

0, Q = Q /3<7a<7/3? 

0, T 2 (g 2 )=3M 2 |/ 2 | 2 , T3(g 2 ) = ^(g 2 -M 2 + m 2 )|/ 2 | 2 + 3g 2 Re/ 1 /*, 
3|g 2 /i + \{q 2 - M 2 + m 2 )/ 2 | 2 , T 5 (g 2 ) = -3g 2 ImA/ 2 *. 



(13) 



B. Expressions for the observables 



Using the definitions of the cross section ([2]), of the leptonic ([3]) and hadronic tensors, 
one can derive the expression for the unpolarized differential cross section in terms of the 
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structure functions W\p (after averaging over the spins of the initial particles) 



da un a 2 p 



dQ 2q 4 W 



(14) 

where r = q 2 / (4M 2 ), t = (fci —pi) 2 and u = (fcj — p 2 ) 2 - In the reaction CMS this expression 
can be written as 



da % 



a 2 p 



dtt 2q A W 



(A + B sin 2 



A = \q 2 h + ~(<? 2 -M 2 + m 2 )f 2 \ 2 , B = 2rp 2 [q 2 \f 1 + f 2 \ 2 - M 2 |/ 2 | 2 ], (15) 

where 9 is the angle between the momenta of the axial-meson (p) and of the electron beam 
(k). Integrating this expression with respect to the axial-meson angular variables one obtains 
the following formula for the total cross section: 



a tot (e + e~ - ira,) = [s^ + \{q 2 - M 2 + m 2 )f 2 \ 2 + 

+4rp 2 [q 2 \f 1 + f 2 \ 2 -M 2 \f 2 {- 



(16) 



Let us define an angular asymmetry, R, with respect to the differential cross section, u n / 2 , 
measured at 6 = tt/2, 

da un 



dn 



a v/2 (l + Rcos 2 6),R= -B/(A + B). 



(17) 



As it was previously shown in the case of e + + e~ — > d + d [17(|, this observable is very 
sensitive to the different underlying assumptions on the axial-meson FFs and does not 
require polarization measurements. 

The differential cross section in terms of FFs fi{q 2 ) contains not only the moduli of these 
FFs but also their interference (1151) . One can choose a linear combinations of these FFs, in 
such a way that the unpolarized differential cross section will contain only moduli. Let us 
introduce new FFs gi(q 2 ) which related to the old ones as follows 



fi= 9i+92, f2 = cg 1 + dg 2 , 



(18) 



where c = —2q 2 / (q 2 — M 2 + m 2 ), d = (q 2 + M 2 — m 2 ) / (M 2 — q 2 + m 2 ). Then the structure 
functions Wi describing the unpolarized part of the hadronic tensor can be written as 

4p 4 q A 



W^q 2 



(q 2 — M 2 — m 2 ) 2 



02 



(19) 



_(g 2 -M 2 + m 2 ) 2 ' y 1 (g 2 -M 2 -m 2 ) 2 
The cross section can be written, in the general the sum of unpolarized and 

polarized terms, corresponding to the different polarization states and polarization directions 

of the incident and scattered particles: 

da da un 

~d~Q = dVt ^ ^ ^^~ >X R^zRzz ~l" PxzRxz ~\~ Pxx(Rxx Ryy) ^-PyzRyz] ■> (21) 

where Pi, Pij, and Rij, i,j = x,y,z are, respectively, the components of the vector, tensor 
polarization and of the quadrupole polarization tensor of the outgoing ai-meson Q^, in 
its rest system and da un /dQ is the unpolarized differential cross section. A is the degree of 
longitudinal polarization of the electron beam. It is explicitly indicated, in order to stress 
that these specific polarization observables are induced by the beam polarization. 

Let us consider the different polarization observables and give their expression in terms 
of the 7* — > din transition FFs. 

• The vector polarization of the outgoing axial-meson, P y , which does not require po- 
larization in the initial state is 



1^ 



(q 2 + M — m 2 ) 2 - AM q sin(20)Im/i/ 2 *, 



(22) 



V 8a 

where a = A + B sin 2 9. One can see that this polarization is determined by non-zero 
phase difference of the complex FFs fi and / 2 . 



• The axial-vector meson can be tensor polarized also in case of unpolarized initial 
beams. 

The components of the tensor polarization are 
3 1 



4 a 



sin^ 9 



q 2 h + \{q 2 -M 2 +m 2 )f^ 



- A - V^1% + M -rn 2 ) 

4 a q 2 I 



q 2 h + \{q 2 -M 2 + m 2 )h 



+ 



+ 



(q 2 + M 2 — m 2 ) 2 - 4M 2 q 2 



_ (? 2_ M 2 +m 2 )|/2|2 + ? 2 Re/i/ * 



}■ 



P — 

J- ZZ 



3 1 1 



8 (j M 2 g 2 



( q ' + M 2 -m 2 ) 2 -4M 2 q 2 



q 2 h + ^q 2 -M 2 + m 2 )f 2 



+ 



+ - sin 2 9 [-g 4 |/i| 2 + ?V + 3M 2 - 3m 2 )ReA/*+ 
+ (2gV-M 2 )-^(g 2 -M 2 + m 2 ) 2 )|/ 2 | 2 ]}. 



(23) 



A possible non-zero phase difference between the axial-meson FFs leads to another 
T-odd polarization observable proportional to the R yz component of the axial-meson 
tensor polarization. It takes the form 



P -J_^L 

VZ ~ 2, 



(q 2 + M — m 2 ) 2 - AM 2 q 2 sin0Im/i/ 2 *. 



(24) 



• Let us consider now the case of a longitudinally polarized electron beam. The other 
two components of the axial-meson vector polarization (P x , P z ) require the initial 
particle polarization and are 

_I S i n {2q 2 {q 2 + M 2 — m 2 )\h\ 2 + \{q 2 - M 2 ) 2 - m 4 l \f 2 \ 2 + 

+ [g 2 (3g 2 - 2M 2 - 2m 2 ) - (M 2 - m 2 ) 2 ] Re/i/ 2 *} , 
11 „ 1 



Px = 



2 ct 



COS0 



q 2 h + - 2 {q 2 -M 2 + m 2 )f 2 



(25) 



III. SPIN DENSITY MATRIX OF AXIAL MESON 



For unstable particles, the vector and tensor polarizations are directly related to the 
angular distribution of their decay products; one can show that the angular distribution 
can be expressed in terms of the spin-density matrix. Let us calculate the elements of the 
spin-density matrix of the axial-meson which is produced in the reaction e + + e~ — > 7r + a±. 
The calculation is done in CMS of this reaction. 

In case of unpolarized initial lepton beams, the convolution of the lepton and hadron 
tensors can be written as 

S un = S^Ut, (26) 

where is the polarization four-vector of the detected axial-meson and the S^ v tensor can 
be represented in the following general form 

Sfu, = Sxg^ + S 2 q»q v + S 3 k lfl k lu + S 4 (k lfl q u + q^k lv ) + iS b (k x ^q v - q^k lv ). (27) 

The functions Si(i — 1 — 5) can be written in terms of the two transition FFs of the axial- 
meson. Their explicit form is 



Si = -q 2 



q 2 fi + 2 (q 2 -M 2 + m 2 )f 2 
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s 9 = 



q 2 {q 2 - M 2 ) - ^{q 2 - M 2 + m 2 + 2Wpcos9f 
+q 2 (q 2 + M 2 -m 2 - 2Wpcos9)Ref 1 f 2 *, 

1 2 



\h? + 



s 3 
s. 



-4 

,4 



q*f 1 + -tf-M 2 + m 2 )f 2 



f 2 + m 2 ) 

f(q 2 -M 2 + m 2 
r2 , m 2\/„ 2 



4 = 2g 4 |/i I + 2g 2 (g 2 - M 2 + m 2 + Wpcos0) 2 Re/i/ 2 * + 

^ .1 1 / - - W - ; i / - Oil / i \ I J 

2| 



+^(g 2 — M 2 + m 2 )(q 2 - M 2 + m 2 + 2Wp cos 



S 5 = -2q 2 Wpcos9Imf 1 f 2 *. (28) 

The T-odd structure function S$ is not zero here since the transition FFs of the axial-meson 
are complex functions. 

The elements of the spin-density matrix of the axial-meson are defined as 



PiijPu 



(29) 



S Pmm > = S, U U^U^\ S = S, u [~g, u + } p 

where S = 2q 2 (A + _Bsin 2 6>), and Ujf 1 ^ is the polarization four-vector of the axial-meson 
with definite (m = 0, ±1) projection on the z axis. In our case it is directed along the 
axial-meson momentum and thus Uj^ are the polarization vectors with definite helicity. 
The elements of the spin-density matrix of the axial-meson are 



P++ 



h {1+cos29) 



q 2 fi + \(q 2 -M 2 + m 2 )f 2 2 



(g 2 + M 2 -m 2 ) 2 |/ 1 | 2 + 2pV + M 2 (g 2 - M 2 + m 2 ) Re/i/ 2 *+ 



+ (m 2 M 2 +p 2 q 2 )\f 2 \ 2 }, 



P+- 



P+o = 



2S 



sin 2 6 



q 2 f l + }-tf-M 2 + m 2 )h 



sin 9 cos 9 { (q 2 + M — m )\fi\ + (1 



q fr . 
+hq 2 - M 2 - m 2 )\f 2 \ 



M 2 - m L 



) (g 2 + M 2 -m 2 )Re/ 1 /*+ 



+ 2p 2 / 2 /r}, 



P-o = -p+o, Po+ = P+o, Po- = P-o- (30) 

The spin-density matrix is normalized as Trp = 1 or p ++ + p + p 0Q = 1. The element p +0 

is complex and the real and imaginary parts are written as: 



Rep+ = -^W^sin#cosfl<j (q 2 + M 2 - m 2 ) \f ± \ 



S V 2 



M 2 



m 



( q 2 -M 2 -m 2 )\f 2 \ 2 + 1 - 



3q 2 - 2M 2 - 2m 2 - 



(M 2 -m 



2\2 



ReA/, 



2 ( ' 



Imp+o = 



2 4 

p f It . 



S 



-sm29lmhr 2 . 



(31) 



11 



Let us consider the case when the electron beam is longitudinally polarized. Then the 
convolution of the spin-dependent part of the lepton and hadron tensors can be written as 

S(X) = S^(\)U„U;, (32) 

where A is the degree of the electron beam polarization and the S^(A) tensor can be written 

as 

S» v (\) = Qie lxva pk la k 2 p + Q 2 (q^a u - q^a^) + Qz(q^a u + q v a^), (33) 

where = e^a^Pakipk^, P — Pi — p 2 - The structure functions Qi(i = 1 — 3) can be written 

in terms of the two transition FFs of the axial-meson as 

1 



-2iX 



q 2 fi + ^(q-M 2 + m 2 )f 2 



q 2 fi + ^q 2 -M 2 + m 2 )f 2 



J 2 i 



Qi ~- 

Q 2 = -2iARe 

Q 3 = 2Ag 2 Im/ 1 / 2 *. (34) 

The T-odd structure function Q 3 is not zero since FFs are complex functions in the TL 
region. 

The elements of the spin-density matrix of the ai-meson that depend on the longitudinal 
polarization of the electron beam can be defined as 

sp mm i\) = s^\)u<r ) u<r' ) *, (35) 

and they are expressed in terms of FFs as 



q 2 fi + ^q 2 -M 2 + m 2 )f 2 



P++( A ) = -P— ( A ) = ^<? 2 cos# 
Poo(A) = p + _(A) = p_ + (A) = 0, 
P + o(A) = -|^ 4 sin^{(g 2 + M 2 -m 2 )|M 2 + 
M 2 -m 2S 



+ 1- 



q 



q^ + M 2 - m 2 )ReA/ 2 * + -(q 2 — M — m 2 )\f 2 \ 



Po+(A) = p%(A), p_ (A) = p +0 (A), Po _(A) = p*_ (A). (36) 

The spin-density matrix element p+o(A) is a complex quantity and its real and imaginary 
parts are 



Rep +0 (A) 



A FT 



q i smel(q 2 + M 2 -m 2 )\f 1 \ 2 + -(l- I ' f ' 



1 

+ 2 
A 



S V 2 

(3g 2 - 2M 2 - 2m 2 - 



q 



(q 2 -M 2 -m 2 )\f 2 \ 2 + 



2 (M 2 - m 2 ^ 2 



Re/:/* , 



Imp +0 (A) = -V2rp 2 g 4 sin0Im/i/ 2 *. 



(37) 
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The axial-meson FFs are complex functions in the TL region. So, for the complete deter- 
mination of FFs it is necessary to measure three quantities: two moduli of FFs and their 
phase difference. Therefore, the measurement of the unpolarized differential cross section 
does not allow to determine completely FFs. It is necessary to determine the spin-density 
matrix elements of the produced axial meson measuring the angular distribution of its decay 
products. 

The measurement of the angular dependence of the unpolarized differential cross section 
allows to determine the structure functions A and B. So, one can determine the following 
ratio 



v 



2M 2 



1 + 2r cos a + | 1 — r z 



M 2 " 
7 



1+1 + 



m 2 - M 2 ' 



if m 

r cos a + - 1 H 



M 2 



(38) 



where r = I/2I/I/1I and a is the relative phase of two complex FFs f\ and f 2 defined as 
follows: a = a>\ — a 2 , where a>\ = arg/i and a 2 = arg/ 2 - Thus, Eq. [38] contains two 
unknown quantities: a and r. Another equation for the determination of these quantities 
can be obtained, for example, from the following ratio of the spin-density matrix elements 
of the produced axial meson 



M. - Sr |^(g 2 + M 2 - m 2 ) 2 + [2p 2 q 2 + M 2 (q 2 — M 2 + m 2 ) 



P+- 

+ (m 2 M 2 +p 2 g 2 )r 2 } 



r cos a+ 



(39) 



1 T 1 
q 4 + q 2 {q 2 — M 2 + m 2 )r cos a + -(q 2 - M 2 + m 2 ) 2 r 2 



And, finally, the phase difference a can be determined by fixing the sign of sin a from the 
measured spin-density matrix element Impo+- Note that the complete determination of 
two complex FFs f\ and f 2 does not require the polarization of the initial beams. The 
measurement of the angular distribution of the unpolarized differential cross section allows 
to determine the structure functions A an B separately. Knowing the ratio r, one can 
determine the moduli and |/ 2 |. 



IV. NUMERICAL RESULTS 



VMD inspired models have proved to be very successful in describing the structure of 
hadrons. Such models contain a small number of parameters, with transparent physical 
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meaning, and can be analytically extended to the full region of momentum transfer squared. 
A monopole-like behavior reproduce quite well the existing experimental data on pion FFs, 

n 

and satisfies pQCD asymptotic [24j. In order to predict the behavior of polarization observ- 
ables, we suggest a simple model for the a\ transition FFs, in TL region. 

We used a simple VMD-based parametrization saturated by vector mesons. The contri- 
bution of one vector meson is given by the Breit-Wigner form 

^- Mrffib? '- 1 - 2 - (40) 

where M v and T v are the mass and the width of a vector meson carrying the interaction. 
In general one should introduce all allowed vector mesons, but, as shown in Refs. 
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251 ]. the largest contribution to the cross section is given by the p(770), p'(1450), and, at 



higher energies, p"(1700). We consider only data for the total cross section at ener gies 



above the a>iii kinematical threshold, which have been compiled from Refs. [5j, [2 



27|, 



. The experimental data 5] show a clear contribution from J/tp around y/s = 3.01 
GeV. We excluded the corresponding (four) data points from the fit. The final form of our 
parametrization of the two 7*7ra! transition form factors 

f C p ,M p C lA M l C^Mf 

J* M 2 -q 2 +iM p T p M 2 -qi + iM p ,T p , M}-q* + iM p »T/ ' ' 1 ] 
Such parametrization has, therefore, in total six parameters: three normalization constants, 
Cp,i, C p ' t i, Cp/'j for each FF. 



The result of the fit for 

7r + 7r~7r + 7r~ process from 



;he available data on the total cross section of the e + e _ — > 
2gI . 27 . 28 ] is shown in Fig. [TJ and corresponds to x 2 l n df = 
455/154 = 2.9. The different sets of data points are quite dispersed, especially in the 
region of the maximum. The resulting normalization constants are given in Table [B The 
7*7ra! transition form factors (moduli) are presented in the FigfSJ Peaks can be seen in 
correspondence with the masses of the chosen vector mesons. For I/2I, one can see a bump 
around q 2 = 4 GeV 2 , which results from the interference of the different terms. 

n 

Such fit assumes that the intermediate state a\T[ saturates the cross section (4|. If other 
intermediate channels contribute to this yield, and only a fraction of the cross section is 
due to the aiii intermediate state, then (assuming no dependence on q 2 ) the normalization 
parameters should be rescaled by the square root of that fraction. 

Figs. [3],|1] show the vector and tensor polarization observables. One should note here that 
the vector polarization of vector mesons can not be measured through their decays which 
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i 








1 


2.35 ±0.01 


-0.089 ± 0.008 


-0.131 ±0.009 


2 


-4.45 ±0.11 


0.56 ±0.05 


0.38 ±0.03 



TABLE I: Optimal values of the normalization constants C P: i,Cp> : i,C p » : i obtained by the fitting 
procedure. 



o 




q 2 [GeV 2 ] 



G. 1: (Color online) Total cross section data for the annihilation e + e — > 7r + 7r tt + tt from 



26, 



27, 



281 ] . The line represents the fit of the total cross section ()16p for the reaction e + + e 



ai + 7r with VMD-based model of transition form factors (jH 



are driven by strong and electromagnetic interaction with conservation of P-parity 29( . Fig. 
|5] shows the predictions for the density matrix elements. These quantities can be quite large 
and show a particular behavior, which can be experimentally verified. Finally, Fig. [H] shows 
the predictions for the ratios Ri, R2 Eqs. f )38|39l) . 



V. AXIAL MESON DECAY TO THREE PIONS 

For completeness, we derive the angular distribution for the decay of the ai-meson to 
three pions 

atip) - 7r+( gi ) + n-(q 2 ) + 7T + (q 3 ), (42) 

where the notation for the particle four-momenta is given in brackets. We assume that this 
decay takes place through a quasi-two-body production of pn and an where a (or / (600)) 



15 




' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' n 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 

q 2 [GeV 2 ] q 2 [GeV 2 ] 



FIG. 2: q 2 dependence of the 7*7rai transition form factors f\, f% (I4ip . 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

e [rad] [rad] [rad] 



FIG. 3: Angular dependence of the vector polarization observables P x , P y , P z at q 2 = 2 GeV 2 , 
Eqs. (I22l25|) . 

is a scalar meson with the following quantum numbers I (J PC ) = + (0 ++ ). The p- and 
cr-mesons are decaying predominantly to tttt final state. 

The matrix element of this decay corresponding to the following mechanism a\ — > (p7r + 
(T7r) — > 3n can be written as 

M = aqi ■ U + bq 2 ■ U = U„A„, (43) 
16 



a 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

[rad] 9 [rad] 



FIG. 4: Angular dependence of the tensor polarization observables P xx , P xz , P zz , P yz at q 2 = 2 
GeV 2 , Eqs. ([23l2iD . 

where U p is the polarization four-vector describing the ai-meson and the functions a and b 
are (in the ax-meson rest system): 

2Mu 1 ) - 2MF 2 u 2 ] + 
wi)] + 

(44) 

Here tui and 1^2 are the energies of the positive and negative pions. The expression for the 
matrix element takes into account the identity of the two positive pions. For the intermediate 
p- and cr-mesons we use a standard Breit-Wigner form and 

D{q 2 ) =q 2 -m 2 p + iT p m p , R(q 2 ) = q 2 - m 2 a + iT a m a , 

where m p (T p ) and m a (T a ) are the masses (widths) of the p- and cx-mesons, respectively. 
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a = G[R-\Q 2 ) - R-\ q *)] + l -D-\q 2 )[{F 2 - F X ){M 2 +m 2 - 

+ ^D-\Q 2 )[(F 2 - Fx)(M 2 - m 2 - 2Muo 1 ) + 2MF 1 u 2 ], 
b = GR-\Q 2 ) - £>" 1 (g 2 )[(F 2 - F 1 )(Mu 1 - m 2 ) + MF X (M - 
+-D-\Q 2 )[{F 2 - F 1 )(m 2 — M 2 + 2Mu 2 ) - 2MF 1 u 1 ), 
q 2 = M 2 + m 2 — 2Mw b Q 2 = m 2 - M 2 + 2M(oo 1 + uj 2 ). 




e [rad] [rad] 9 [rad] 



FIG. 5: Angular dependence of the elements of the spin density matrix poo, p^ , p+o at g 2 = 2 

GeV 2 , in case of unpolarized collisison (solid lines) from Eqs. (|30|3ip . In case of longitudinally 
polarized electrons the corresponding elements are shown also (dashed lines, (d),(e)) from Eqs. 
([37D and (f) from Eq. ([36}. 




.-3 I i i i I i i i I i li i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i 

1 2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 



q 2 [GeV 2 ] q 2 [GeV 2 ] 

FIG. 6: q 2 behavior of the ratios R 1: R 2 , Eqs. (I38T39D . 
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The quantities G and F± i2 are defined as follows 

G 9aair9cmiT i F*i 9pirnfii i 1) 2, 

where f\, i = 1, 2 are the coupling constants defining the a\ — > p7r vertex. 

The definition of the different coupling constants and their relation to the corresponding 
decay width is given below. 

The a — » irir decay. The amplitude of this decay can be written as: M = g^^l^fp, 
where <pi and ip are the wave functions of the pions and cr-meson, respectively. Then the 
expression of the decay width in the a-meson rest frame is 

r( ( T^7T7r) = -|^ v /m2-4m 2 . (45) 
167rm^ v 

The p — > 7T7T decay. The matrix element describing this decay can be written as: M = 
Qp-K-K^pisii — q2)^flf 2 , where e M is the p-meson polarization four-vector. qi(q 2 ) and (pi((p 2 ) 
are the four-momentum and the wave function of the first (second) pion respectively. The 
expression for the width of this decay in the p-meson rest frame is 

r(p--) = J^K-4m 2 ) 3/2 . (46) 

The a\ — » an decay. The matrix element in this case is: M = gw^u(?i — 12)^*^*, where 
11^ is the ai-meson polarization four-vector and qi(q 2 ) is the a(ix— )-meson four-momentum. 
The width of this decay in the ai-meson rest frame is 

r(oi -> an) = /%-[(M 2 + m 2 - m 2 ) 2 - Am 2 M 2 f 2 . (47) 

The ai — > pn decay. The matrix element of this decay is determined by two coupling 
constants and can be written as: M = fi(q 2 U ■ e* — q ■ Uq ■ e*) + f 2 (q-p 2 U ■ e* — q ■ Up 2 ■ e*), 
where e^{U^) is the polarization four-vector of the p(ai— )-meson and q(p 2 ) is the four- 
momentum of the p{n— )-meson. The expression for the width of this decay in the a\- 
meson rest frame is 

2 

77? 

Tia^pn) = y^-^[(M 2 + m 2 - m 2 ) 2 - 4m 2 M 2 ] 1/2 (ci|/i| 2 + c 2 \f 2 \ 2 + 2c 12 Ref 1 f 2 *), 
ci = 8m 2 p M 2 + (M 2 + m 2 - m 2 ) 2 , 

c 2 = 3M 2 (2m 2 - M 2 ) + (m 2 - m 2 ) 2 + 2^-(M 2 - m 2 ) 2 , 

p m 2 p 

Cl2 = M 4 + 4M 2 (M 2 - m 2 p - m 2 ) - (m 2 p - m 2 ) 2 . (48) 
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Let us calculate the square of the matrix element for the a\ — > 37r decay using the helicity 
formalism. In the general case the ai-meson is described by the spin-density matrix and 
one can write: 

\M\ 2 = p xx ,B xy , B XX , = A X A* X „ (49) 

where A x = U^A^ and Uj^ is the ai-meson wave function with definite helicity, p X y is the 
ai-meson spin-density matrix in the helicity representation. The standard choice of the z 
axis is along the normal to the decay plane in the ai-meson rest frame (we chose it along 
the direction q 2 x <fi, where 91(92) is the three-momentum of the first (second) pion). Since 
all pion momenta are perpendicular to the z axis the quantities B xx > with zero helicity A or 
A' are equal to zero. Finally, one obtains: 

|M| 2 = p++B ++ + p__fi__ + 2Rep+„Re J B+_ - 2Imp + _Im J B + _. (50) 

As it is known, the decay of a particle into a three-body final state is characterized by two 
independent variables. We choose the energy of the negative pion u 2 and angle 6 between 
the momenta of the first and second pions, i.e., between 9! and 9*2. The x axis is directed 
along q 2 . In this case the quantities B xx > are 

B±± = ^\a\ 2 + q^\b\ 2 + 2 qi q 2 (cos9Reab* T sin 9lmab*)}, 
ReB + _ = -i[gt|a| 2 cos2# + g||&| 2 + 29192 cos0Rea&*], 

Imi? + _ = — 9192 sin^Reafe* — -9^|a| 2 sin20, (51) 

where 91(92) is the value of the three-momentum 9*1(92). 

The energy and angle distribution of the decaying ai-meson is described by the following 
expression 

dr(ai^37r) 1 9i92 r?l/f . 92 /)i-iu#|2 , KO \ 

= - — — \M — too H kjicosft \M\ . 52 

duj 2 d6 (2tt) 4 AM 1 91 J 

The energy u\ is not an independent variable and it can be determined in terms of two 
independent variables using the following identity 

29! 92 cos# = M 2 + 2m 2 - 2Mu 2 + 2wi(cj 2 - M). 

We introduced above the general spin-density matrix for the description of the polariza- 
tion state of the ai-meson. In the coordinate representation its expression is given by Eq. 
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In the ai-meson rest frame this formula is written as 

1 i 

pij = ~^£ijkSk + 3Qij, ij = x, y, z. (53) 

This spin-density matrix can be written in the helicity representation using the following 
relation 

Pay = Pij ul X) *uf\ A, A' = +, 0, (54) 

where [// are the ai-meson spin functions which have the ai-meson spin projection A on 
the quantization axis (z-axis). They are 

tf (±) = =F^(l,±i,0), U^ = (0,0,1). (55) 

The elements of the spin-density matrix in the helicity representation are related to the ones 
in the coordinate representation by such a way 

113 1 3 

P±± — g ± — -jQzz-, Poo — g + 3Q ZZ , p + _ = — -(Qxx — Q yy ) + 3iQ xy , (56) 

13 13 

P+o = ^^{s x -iSy)--j={Q xz -iQ yz ), p_o = ^/|(s a; +^j / ) + -^=((5 xz +2(5j /z ), Paa' = (pa'a)*- 

To obtain these relations the condition Q xx + <5 ra + Q zz = was applied. 

The factor 1/3 in the unpolarized part of the spin-density matrix was introduced, as for 
the decay it is necessary to average over the spins of the ai-meson. 



VI. CONCLUSION 

Using the parametrization of the electromagnetic transition 7* — > ai(1260)7r in terms of 
two FFs, we investigated the polarization phenomena in the annihilation reaction (JTJ). We 
calculated the differential (and total) cross section and various polarization observables as 
functions of the FFs. The spin-density matrix elements of the produced axial meson have 
been also calculated. Explicit formulae for the decay of the ai-meson into three pions have 
been given. 

FFs are complex in the time-like region and have been parametrized according to a VMD 
inspired Q 2 dependence, saturated by vector mesons. The parameters have been fitted to 
the data, assuming that all the four pion yield is due to reaction ([T]). If other intermediate 
channels contribute to this yield, and only a fraction of the cross section is due to the aiir 
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intermediate state, then the normalization parameters should be rescaled by the square root 
of that fraction, assuming no dependence on q 2 . 

The q 2 dependence of the FFs shows, as expected, peaks in correspondence of the masses 
of the considered vector mesons. An interesting feature appears for f 2 at Q 2 ~ 4 GeV 2 , 
which is due to an interference of the different terms. Note that none of the considered 
vector mesons corresponds to a mass of such value. An interference of such origin can be 
the reason for the plateau seen in the cross section data, which, however does not appear 
with large intensity in our fit. 

The reaction (TjQ) has been clearly detected in the experiments. The present results can 
be useful for the analysis of the experimental data and for the determination of the 7* — > 
ai(1260)7r transition FFs. 
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